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MERCURY'S HELIUM EXOSPHERE AETER MARINER lO's 
THIRD ENCOUNTER 

S. A. Curtis and R. E, Hartle 
Laboratory for Planetary Atmospheres 

ABSTRACT 

From a comparison of the Mariner 10 third encounter UV 
spectrometer data with intensities generated from a newly 
constructed model exosphere, we have derived a new value 
of 4, 5 X 10 for the fraction of the solar wind He'*"*' flux 
to be intercepted and captured by Mercury's magnetosphere 
if the observed He atmosphere is maintained by the solar 
wind. If an internal source for He prevails, the corres- 
ponding upper bound for the global outgassing rate is esti- 
mated to be 4, 5 X 1Q22 sec-l. These values differ from 
those given earlier due to the present use of a surface tem- 
perature distribution satisfying the heat equation over Mer- 
cury's entire surface which employs Mariner 10 determined 
mean surface thermal characteristics. We also use the 
mean stand off distance of Mercury's magnetopause averaged 
over Mercury's orbit. We find the agreement between the 
observed and calculated intensities to be good and believe that 
the minor discrepancies that exist on the nightside of the 
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ABSTIUCT 


lormimitor are explicable in terns of differences between 
actual and computed surface temperatures and the resulting 


scale hetg'lit structure clianges. We attribute tliose tem- 
perature differences to the inhomogcnelty of the pliyslcal 
properties of the surface of Mercury, 
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MERCURY'S HELIUM EXOSPHERE AFTER MARINER lO's 
THIRD ENCOUNTER 

INTRODUCTION 

As a result of Mariner 10' s three encounters with Mercury, the prevailing 
view of the relatlonsMp between the atmosphere and the solar wind has changed 
dramatically. A review of the observations and Interpretation of Mercury's at- 
mosphere has been presented by Kumar (1976), Our Interest here lies primarily 
In the Interpretation of the third encounter observations of the He atmosphere at 
Mercury's evening terminator (Broadfoot et al. , 1976a), We will compare our 
theoretical 584 A Intensities to the Intensities observed by Mariner 10 In order 
to ascertain the validity of the model's assumptions concerning the He atoms' 
accommalatlon with local surface conditions. Our investigation is constrained 
by uncertainties in the derived surface temperature due to inhoniogeneltles In the 
surface properties. We will show that the data appear to support the model's 
assumptions to a degree wMch allows new updated estimates of the total fraction 
of solar wind He++ captured by the Hermian magnetosphere and which impacts 
the planet's surface, maintaining the observed He exosphere solely by this exter- 
nal source. Also, a new upper bound to the total surface outgassing rate will be 
given. 

Pre-encounter models (Hartle, et al. , 1973; Hodges, 1974) assume a direct 
interaction of the ambient solar wind with the Hermian surface, m those models 
the solar wind provided the source for the atmosphere, and its magnitude was 
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detei'mlnud by the solar wind flux Intercepted by Mercury, The Mariner 10 en- 
counters discovered tliat direct solar wind Interaction was generally not possible^ 
since Mercury was found to possess a magnetosphere (Ness et al. , 1974} Ogllvle 
et al . , 1074). This magnetosphere deflects most of the solar wind from Mercury 
under typical solar wind conditions (Slscoc and Christopher, 1975). The encount- 
er measurements also revealed a He atmosphere (Broadfoot et al. , 1974). From 
encouutor measurements Hartle et al, (1975) construeted a model to account for 
the observed magnetospherlc shielding of Mercury's surface from the solar wind. 
The model Indicated the observations Implied a tenuous He atmosphere which is 
colUslonless down to the planet's surface and lienee Is an exosphere. Fx’om the 
model, preliminary estimates were made for the magnetospherlc capture fraction 
of solar wind He'*“*‘ required to maintain the observed He exosphere. In addition, 
an upper bound for the surface outgasslng rate required to maintain the He exo- 
sphere was also obtained from the model. For simplicity, the early model as- 
sunied a surface temperature distribution based on radiative equilibrium. This 
assiunption led to a good approximation of the temperatures over most of the day- 
side, but was not appropriate near the evening and morning terminators (Hartle 
et al. , 1975). The need for a more realistic temperature distribution Is evident 
from the Mariner 10 third encounter UV observations of the sensitive terminator 
x’cglon (Broadfoot et al, , 1970a). A minimum requirement for a surface tempera- 
tux'e distribution in the present case would be for it to satisfy the heat equation 
over the entire Hexunian surface using the surface thermal properties determined 
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by Mariner 10 (Chaae et al. , 1970). The differences in global parameters such 
as total atmospheric content are not largo when the temperatures are determined 
from either the radiative oqulliblrum condition or tlie heat conduction equation. 
However, the localized UV measurements of Mariner 10 are extremely sensitive 
to tlie terminator temperature distribution and hence a very realistic temperature 
model is needed if theory is to bo quantitatively compared with experiment in this 
region (liar tie et al. , 1975), 

METHOD & CALCULATION 

In this section we discuss in detail the three stages of calculations that are 
necessary in order to generate a model helium exosphere for Mercury from which 
detailed comparisons with Mariner 10 encounter UV observations are possible. 

The first of these stages is the calculation of the surface temperature distribution, 
obtained by solving the heat conduction equation over the planet’s surface. Tlie 
second stage incorporates the surface temperature distribution into a Monte Carlo 
model for the surface helium exobase density normalized by a constant factor. 

The tliird and final jstage involves a solution of the collisionless Boltzmann equa- 
tion for the global He density distribution in the exosphere, where the highly non- 

J: 

uniform exobase temperatures and densities provide the needed boundary conditions. 
We now elaborate on each of the three parts of tlie model generation. 

In calculating the exobase temperature, T(z, t), we solve the heat conduction 

equation (Bichtmyer, 1962), assiuning a homogeneous surface 

3T(z.t) = K ^(z.t) " (I) 

~ir~ 3z2 
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where p 1b the surface density, e'ls the specific heat capacity, K is the thermal 
conductivity, t is the time and z is tlio distance below Mercury's surface. In 
solving (L) we applied tlio following boundary conditions! 

K DT(z,t) 

K aT(z.t) 

3z 

whore e is the emlsslvlty, a is tho Steplian-Boltznmnn constant, a is the albedo, 
and Zq is the depth of the bottom of the slab In wliich the solution Is sought. P(t) 
is tho solar Insolation as determined from Mercury's ephemeris. In (L) and (3) 
we have neglected the heating contribution due to radioactive decay of uranium 
and thorium In the Hermlan Interior as it is negligible In the case of a luilform 
source of decay under the planet's surface. If, however, the sources of radio- 
active decay are highly non-unifoun, then local hot spots may be produced on the 
night side. The deviations between such hot spots and the luiiform temperature 
model would be similar to those deviations caused by higher thermal Inertias than 
the uniform value assumed. In Table 1 we give the values of the above parameters 
used in the solution. The values characterizing the Hermian surface are taken 
from those adopted by the IE team from their rachometer observations (Chase et 
al. , 1976), The details of the numerical analysis of (1), subject to (2) and (3) as 
well as the analytical form of P(t) have been discussed extensively (Llnslcy, 1965; 
Ingraoetal., 1965; Linsky, 1966; Morrison, L969; Morrison, 1970), In solving 


= 0 




(3) 


= caT4(/=0,t)-(l-a)r(0 


z=0 


( 2 ) 
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(L) wo Imvo choson a grid with a time Inorement At I earth day and a corre- 
sponding depth Increment Ak «« 3 cm as necessitated by stabllUy requirements. 
We chose as a first guess to the solution the radiative equilibrium solution wo 
employed earlier (Hurtle et al, , 1976), Wo note however that particular care '' 
must be taken In numerically solving (L) by the methods cited. Speclflcallyj the 
solution Is prone to oscillatory problems that can destroy the convorgenoe In the 
low tenilperature gradient regions on the night side and In areas With very low 
solar Insolation near the terminators at higher latitudes. The tendency Is to over- 
shoot the true solution. To <;''lmlnate this problem wo apply a dampening factor 
to the computed change of tlie surface temperatures at each Iteration step In these 
oscillation prone regions. We then obtain a smooth convergence with an estimated 
precision of approximately 1°K, The iteration Is performed over multiples of two 
Hermlan years (two orbital periods), the time over which a surface point expnx’l- 
ences a complete solar Insolation cycle. The Insolation received by the aphelion 
Bubsolar point and its resulting temperature response Is shown In Figure 1 for 
the complete two Hermlan year (s«i70 earth days) solar Insolation cycle. The 
unusual form of the Insolation Is derived from the 3/2 spln-orblt resonance of 
Mercury and its relatively high orbital eccentricity. Although wo have derived 
the exobase temperature distribution over the entire Hermlan solar Insolation 
cycle, our interest here Is restricted to the aphelion temperatures at the time of 
the Mariner 10 encounters. In Figure 2 we give the aphelion exobase surface 
temperatures T{\,<f>) as functions of latitude \, and longitude (p , with ip~0° 
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reproBcntIng the noon meCviiian und X ~0® the oquatoi*. In Figures 2b-c \ve show 
eross sections ol this contour plot at low and high latitudes with Uio horizontal 
axis scaled by co^ In order to t^ive the reader a feeling for the longitudinal tem- 
perature gradients (circumference cosX ). As Is readily apparent from these 
two figures, tliere Is a marked asymmetry between die morning und evening ter- 
minators with the smaller evening tempci-ature gradients corresponding to night 
time cooling and the morning temperature gradients being quite steep duo to the 
step function- like turn on of the solar Insolation at sunrise. Also, for large 
distances about the subsolar point on the dayslde, die solution Is very nearly that 
given by radiative equlllbrlimi with T( X , ) -’ 575 f SlnX eos</»J The surface 

temperature on Mercury has perhaps the largest contrast In die solar system. 

At aphelion It spans temperatures from a nlghtslde pole value of to a sub- 

solar value of 

We use the aphOllon exobaso temperature described above to generate the 
exobase helium density distribution. The method we employ Is a Monte Carlo 
technique wherein single He atoms are followed over the planet's surface until 
they either escape thermally or, as Is less likely, are Ionized and picked up by 
the surrounding magnetosphere, The details of the calculation have been presented 
in an earlier paper (Hartle et al. , 1975). The resulting aphelion exobase He den- 
sity distribution n(X , 0, z~0), normalized by the subsolar point density no , is 
shown In Figure 3a. In Figures 3b-c we display cross sections of Figure 3a at 
low and high latitudes, with the horizontal (0) axis scaled by cosX In order to 
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iir order to accurately show the longitudinal surface density gradients. Tlio tig- 
nature of the exobase density response to the exobase temperature distribution 
is apparent from the Steeper gradients In the morning than in the evening. We 
note, however, that non-local effects of a global nature are also apparent from 
an inspection of Figures 3a-e, The maximum density bulge occurs off the etiua- 
tor and marks the most effective trapping region for helium on tlie nightstde, 

The effectiveness of the trapping is determined both by how cokl it is and the ex- 
tent of tlio cold region on the surface. These effects combine to determine the 
length of time an atom will spend balltstically hopping around In a given region 
and hence gives a measure of how long it remains trapped. Thus the maximuni 
density ^ o-.'.x vot appear near the poles which are tlie coldest areas on the planet 
nor does it occur at the equator which has the largest low temijerature region. 
Rather, the maximum density occurs in between, where both of those characteris- 
tic trapping parameters combine to form the best trap. The smallness of the 
polar trapping region is readily apparent from both Figures 2a-c and 3a-c. Most 
of the nlghtside structure seen in Figure 3a can be attributed to global processes 
and not the relatively featureless local nlghtside exobase temperatures of Figure 
2 a. 

Using the exobase temperatures and He exobase densities described above, 
our model liermian exosphere is generated directly as a solution to the collision- 
less Boltzmann equation with these exobase distributions as boundary conditions. 
The details of the method of solution of tills equation for nonuniforni planetary 



exobaec conditions is given by lUu’tlo (197L), liatlier than computing the entire 
exosphere we iiave, I'or reasons of economy, rcGtrlcted our calculations to those 
necessary to compute tlio coliunn densities needed for comparison with the Mari- 

\i 

nor XO tlilrd encounter UV observations. The column derisltles wlUeh wo compute 
are of tiie form 



where So is the intersection of the viewing direction with the planet when sucli an 
intersection exists (for views across the limb Sq = and Sc Is tlio spacecraft 
position at the tlnie of viewing. The instriunent viewing slit (Broadfoot et al. , 
I97Gb) is divided into segments and for each segment a coliunn density of the form 
of (4) is computed. The integrated contribution of these segments for a given 
slit (accounilng for the instrument's transmission function) gives tlie normalized 
intensity seen by each silt. In order to determine the normalization factor for 
our ■'omputed intehsitles, they are renormalized to those measured off the day- 
side limb by Broadfoot et al. (1976a). In I^igure 4 we compare the resulting re- 
normalized model intensities with Ui_e Mariner iO measured He 584 A intensiies. 
The plot is conveniently broken Into three parts. These parts are: (1) the region 
off tlie dayside limb wiaicli sliows very good agreement with tlie model, (2) the i 
region between the dayside limb and the terminator which shows the observed .^n'f 
tensities exceeding the model by a large amount due to surface scattering of 584^ 
radiation off .the planet's surface, and (3) tlie region to the niglitside of the termin- 
ator whioli follows the data but shows some discrepancies. 
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DISCUSSION n , i 

Returning to Figure 4, wo note the agreement with the observutlcns is good 
off the daysidc limb (region 1) as is expected since radiative equilibrium holds 
on tlie. surface there and thus the local temperatures are nearly Independent of 

( i , 

the local thermal properties of the soil, in contrast, there are notable differen- 
ces on the nlghtside of the terminator (region 3), The large differences between 
model and data in the Intermediate region 2 are attributed to Mercury's UV sur- 
face albedo. We find these night side (region 3) discrepancies to be much less 
than those calculated from a model exosphere by Broadfoot et al. (197Ga), where- 

<' I- 

in a large excess of density was derived in this region. From Figure 4 the nature 
of these differences in region 3 are seen to be like a phase shift of the model in- 
tensities to later local times with respect to the observed Intensities. We suggest 
the nlghtside differences near the terminator are due to the extreme sensitivity 
of the intensities to the surface temperatures in this region of large longitudinal 
temperatur^e gradients. The temperatures themselves are sensitive to the local 
thermal parameters and hence will reflect soil inhomogenities. This sensitivity 
arises from the location of the exospheric densities contributing ir ost to the ob- 
served Intensities. Such densities are generally located at several scale heights 

above the surface of Mercury, where Rj is the Ilermian radius and 

GM^ m - 

“ / kT(X,(?5)’ 


(5) 


( 6 ) 




“T 


T 




r 
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1b u ineuMure of tiu.* Ik* utomH ratio of Kt uvltutionul poli'ntial energy to kinetic 
energy in temiB of the gravitational eonBtant G, the maHs of Mercury Mj , the 
lie atom'B niasB ni ami Holtzniann'B con»tant k. The terminator ivgion is charac- 
terized by II '<K)-100 km. Surface clenHiticB do not contribute to the observed 
intensities here due to the shadow of Mereur)' that blocks light to thissuiiace re- 
gion. The densities in this night region arc well approximated by the barometric 
law 

n(X.0) « n(X,0, z»0) exp (-jj— | (7) 

7 

We note that if there exists an error in the exobase temjxjrature, not only is 

the base density distribution modified, but the height distribution is nonlinearly 

o 

altered with a being replaced by a' « 7TTI . Thus, if the actual temperature at 
the exol>ase is lower than that calculated from the heat equation, or' will increase, 
II will decrease and the density will decline more rapidly with height, resulting 
in lower column densities and intensities. To lx.* more quantitative, consider the 
fractional change in density at height z causeil by a difference Ixtween model tem- 
|x*ratures, T,,, , anJ actual temperatums of §T * T>\ - Tmt 


An 

n 


II 




exp 


,a/. 61' 
R. T 


I -1 


(S) 


where we have used the fact that the exobasc densities vary with the exobase tem- 
peratures roughly as T"^ . This density variation with temperature is deduced 
from the minimization of 
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where the Integration is over tlie Hermlan suri'ace. We solve the equation Itera- 
tively for X such that 5 (x)« I and find x - 3. 1(36 ± , 002, For purposes of Illus- 
tration here we will simply use k ~ 3. Now we consider the case In which 5T<0. 

In particular, from Chase , et al. (1976) in their Figure 6 we note that there Is a 
large deviation between model temperatures and observations near 1800 Mercury 
local time; l. e. , near the evening terminator. Specifically, with a thermal in- 
ertia of 0.0020 cal cm-2 sec"l/2 ®K-i similar to that adopted in our temperature 

|6T 

calculations and with the Same type of temperature model, deviations, 0.2 

5T 5 n 

appear to be implied. Using the value ~^ = “0.2 the fractional deviation — at 

2 scale heights above the 'surface is .45, Thus the model would appear to 

overestimate the density by almost a factor of two. We conclude from these re- 
marks that the non-linear amplification of model-actual temperature difference 
effects via the scale height dependence of the 'column density contributions can 
explain the observed phase shift-like differences In the terminator region. We 
attribute the temperature differences to the idealization made in regarding the 
Hermian surface as homogeneous and uniform in terms of its thermal parameters. 
The sensitivity of the surface temperature to these assumptions is readily seen 
from Figures 6 and 7 of Chase et al. (1976). We emphasize here that the measure- 
ments made In the terminator region are Indirect indications of the magnitude of 
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the largo nlghtsldo densltloa, similar to tliose measured directly on the Moon by 
means of a surface mass si>cctrometer (Uodges and Hoffman, 1974). It would 
apixiur that by such direct In situ measurements of density a separation of nlglit- 
slde temperature and densities could be much more precisely obtained. A sAJullar, 
although less precise, separation could be made If the temperatures were mea- 
sured globally by an orbiting satellite allowing the construction of an empirical 
' temperature model reflecting the surface Inhomogeneltles. This motlel coupled 
with UV Intensity observations would separate densities and temperatures better 
than the current encounter measurements of Mariner 10. The differences between 
the homogeneous model Intensities and the Mariner tO observations have been 
taken by Broadfoot et al. (L07t5a) to Indicate departures from the surface physics 
assumed by these models (Hodges, L974; Hartle et al. , 1975). I'rom the results 
obtained here wo feel that for the specific case of helium on Mercury sucli devia- 
tions could well be small since the deviations from homogenlety of the surface 

thermal properties as determined from the Mariner 10 IK experiment (Chase et / ~ 

// 

// 

al. , 1970) appear to account for the differences between model and W data inten- 
sities. A determination of any non-negllglble surface physics deviations would 
require more detailed measurements as discussed earlier. ; , 

From the preceding comparison of the Mariner 10 UV results for atomic 
helliun and oui 'model, we conclude that the agreement Is sMliclently good so as 
to allow a calculation of the total atmospheric content. We are then able to esti- 
mate the capture fraction of solar wind He+''' of the llermlan magnetosphere needed 
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.to maintain the observed atmoaphore. In tiio opposite limit, we also compute an 

I 

i; 

upixjr bound for the HormUm outgassing rate presumably for an Internal radio- 
genic source. In calculating the total atmospheric content, N -p , we integrate 
over a semi-infinite shell containing the Hermlan surface 






/ * r+1 r 2 n 

t i-i““ J 0“^ 


n(u,?i,z) 


whore u = cos X , We find to be well approximated by 

* dll r H(u,^), n(u,0,z=O) 

Jo 

where H(u, 0 ) is the local scale height. The resulting total exospheric content 
we compute Is N'p - G x 10^^ He atoms. This value of N is less than that deter- 
mined earlier (Hartle et al, , 1975) due to the eniianced post sunset evening tem- 
peratures of the ^leat conduction equation solution as sliown in Figure 2, These 
are to be contrasted with the pre-dawn morning temperatures which closely follow 
the radiative equilibrium solution. The resulting evening density depletion of 
Figure 3 reduces the total nightside density by ~,3, As the night side density 
bulge is almost the entire content of the exosphere, the total atmosphere content 
is reduced by a corresponding amount. Also from a comparison of model intensi- 
ties to observed intensities we derive a subsolar point density of Oq = 1 x 10^ 
atoms cm , In obtaining n^ we have used the relation 

in terms of the currently accepted value of g-factor (Broadfoot et al, , l97Ga), 
where I is the observed Intensity and r, is the computed model column density 
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i'rom uqiiifUlon (4). This value of Hq Is In agreement with our previously reported 
value (Hurtle ot al. , 1975) when the differences In observed Intensities between 
Mariner 10 encoiuiter I and encounter III are accounted for. This Is to be expec- 
ted since ny Is determined off the dayslde limb whose surface temperature Is 
nearly indeijendenfc of soil thermal properties. Our present value of ny does 
differ substantially from that calculated by Broadfoot et al. (197(3a) of 4. 5 x ICP 

- I'/Vi 

enr^ , Comparing their Figure 2 to our Figure 4 we note that off the dayslae 
limb the Broadfoot model Is consistently above the data whereas our model trades 
the data quite well in this region. If the Broadfoot et al. model Is normalized to 
ours just off the limb, their value of ny becomes 2. 5 x 10^ enr^ , In agreement 
with the value presented here. VVe now proceed to calculate the capture fraction, 
fy, required for the lie exosphere's maintenance. Slscoe and Christopher (1975) 
have discussed the dependence of the solar wind stand off distance, r^t . at Mer- 
cury as determined from the balance of solar wind pressure and the Hermian mag- 
netic dipole field pressure. From their values of rjjt at aphelion and perihelion, 
the dependence of the dipole magnetic field pressure, ou J-'st ^nd the de- 

pendence of the solar wind density on solar distance R via the continuity equation 
for solar wind He+'^ , we find that the capture cross section, Osc , Varies as 

Ogy - 6.657T R| 6 

where 5 is a geometrical factor which is the squared ratio of the magnetopause 
distance at 90® to the standoff distance 5 = (rst(90“)/(rsi(0°))2 and R is in AU. From 
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FlgLiro 2 of Noss et al, (t97(3) ono obtains un estimate of 6 = 1. 60. Now the moan 
rosltlence tbiie r of an Ho atom on Meroury, before it is lost by an escape process 
is ~ 1,3 X see slightly less than the value reported earlier (Hartle et al. , 
1075). TlUs is on the order of the temperature cycle period of 2 Hermlan years. 
We thus expect the total He exospheric content to be approximately constant with 
variations in the exobase temperatures affecting the shape of tlie exobase density 
distribution. We can then write 

!i 

N-r - 6.65ir iioVtj < R2 (L2) 

fi 

where noUo= 1.5 x 10^ He’^ cm~2 sec-1 is the solar wind flux at 1 AU. The 
n-4/3 dependence reflects both the radical dependence of the solar wind flux on 
distance (R"2 ) and that of the captux’e cross section (K2/3), The bracket < >orbit 
represents an average over the Hermian orbit, 

/ 2it 

(Id [1 +.206 cost?] -2/3 (1-3) 

U 

From (13) and (12) we obtain fc ^ 4. 5 x lO'^l as the Hermian magnetosphere's 
capture fraction. The upper bound on the total outgassing rate^ .^^T, is given by 

T 

(12) directly and is ~4.5 x lo22 He atoms sec-1 , 

This upper bound on the total outgassing rate may also be viewed as tlie upper 
bound to thermal loss. The value reported here is substantially less than that 
given by Broadfoot et al. (1976a) of 3. x 10 23 sec -1 , We attribute this diffei'ence 
to large differences in exospheric scale height structure betxveen the models. 


15 



This can bo seen in the large dlffei’oncos between tlio model Intensities from day- 
side limb to evening terminator in Broadfoot et nl. Figure 2 and our Figure 4, 

In particular the Broadfoot et al. model sliows a much larger bulge in the sunlit 
region near the terminator. We bollevo this bulge is unrealistic as it contains 
model values greater than observed data points which possess surface albedo 
contributions. 

SUMMARY 

From a solution of the heat conduction equation over Mercury's surface we 
have calculated the surface exobase temperatures at aphelion. Wltli tiiese tem- 
peratures we have generated normalized e.xo)jase densities on the surface via 
Monte Carlo techniques and computed the normalization factor by comparing 
model 584 A Intensities with those intensities observed by Mariner 10 at its third 
Mercury encounter. Using these exobase temperatures and densities we have 
computed the total atmospheric content and from this value estimated the total 
magnetos pheric capture fraction of solar wind He+''' needed to maintain the obser- 
ved exosphere solely by this external source, We foujid fg'" 4. 5 x 10“'* . An 
upper bound to the surface outgassing rate of 4. 5 x iO^^sec"* was also derived. 
These results represent a synthesis of the data obtained by the UV, IR and mag- 
netic field Instrimients during Mariner iO’s encounters with Mercury. In closing 
we note once more the sensitivity of the observations to tlie surface temperature 
distribution and the corresponding nonlinear amplification of temperature differ- 
ences in computed model colui/ih densities on the nlghtside near the terminator. 
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Those tompeniUu’o dlCferoncos In tui*n aro doponclent on the Inhomogonlotlcs oi‘ the 
aurfaco properties away from the vicinity of the subsolar iwlnt where radiative 
equilibrium holds, A more refined analysis will most likely require the informa- 
tion that can only bo provided by an orbiting spacecraft which would allow detailed 
observations of the surface temperature at all local times and a wide range of la- 
titudes as well as more detailed density measurements or preferably a lander 
capable of in situ measurements. This contrasts to the three thin slices of en- 
counter data that Mariner 10 has provided us in our first close-up view of Mer- 
cuiy. 
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Table I 


I lCal Comluclioii liquation Parameters 


i’araineter | 

Value 

p (surrace density) 

1.50 gcnr^ 

c (spccirie heat capacity) 

0.20 cal deg-J g'l 

K Cthermal conductivity) 

1.30 X 10'^ cal cm** .sec** deg** 

e (emissivlty) 

0.90 

a (albedo) 

0.10 

a (Stcphen-boltzinann constant) 

1.35 X 10'*^ cal cnr^ deg*'* sec* 
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Temperature (*K) SoUr Inso^atnn r (t) cal cm ^*ec ‘ 


(A) 



Elar%ad Time From Aphelion (Earth Days) 


Figure I (u) Solar insolation history ol the aphelion suhsolar point of Mercury as a function 
of elapsed time from aphelion in earth days. (b)Temperaturehistory of the aphelion sub- 
solar point in Mercury as a function of elapsed tim»* from aphelion in earth Jays. 
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I'igiire 2 (a) llennian siirfaa' temperature disirihution at aphelion asa function of latitude, 
X. and longitude, <j>. (b) Longitudinal cross section of Figure 2(a) for polar latitude, X='’7. 5°. 
^c) Longitudinal cross section of Figure 2 (a) for equatorial latitude, X=2.5“. 
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Figure 3 (a) Hennian surface He density at aphelion as a function of latitude, X, 
and longitude, 0. (b) Longitudinal cross section of Figure 3 (a) for polar latitude, 
X=77.5°. (c) Longitudinal cross section of Figure 3 (a) for equatorial latitudes, X=2.5®. 



Terminator 


Limb 


f'igurc 4 ('omparison of the He 5X4 A evening terminator intensities observed by 
Mariner 10 at its third encounter and model Me exosphere computed intensities. 





